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Abstract 

By using canonical Monte Carlo simulation, the liquid-vapor phase diagram, surface tension, inter- 
face width, and pressure for the Mie(n,m) model fluids are calculated for six pairs of parameters 
m and n. It is shown that after certain re-scaling of fluid density the corresponding states rule can 
be applied for the calculations of the thermodynamic properties of the Mie model fluids, and for 
some real substances. 
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Among the intermolecular effective interaction potentials, the Lennard- Jones one is by far 
the most widely used for approximatin g th e physics of simple nonpolar molecules in all 
phases of matter [l[ [2, [], 4], 5, [f3, Q, |g[ [qL lloL 11, [l2]. The attractiveness of the LJ model is 
mainly due to its more convenient mathematical form than to its accuracy in representing 
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or 



the properties of real fluids. Some modifications of the LJ potential, like exp-6 
the family of Mie(n,m) potentials [13I . Q, have shown to be useful for the description 
of thermodynamic and dynamic properties of some real substances. The Mie(n,m) pair 
potential, which is just a general form of the LJ model, is defined as, 
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where r is the interparticle distance reduced by the particle diameter, which is chosen to 
be the unit length, a = 1; e is the well depth. The temperature of the system is defined as 
T = k B T/e. 



Recently, both theory and molecular simulations have been used 



10 



11 



12 



15 



16 



17 



18|, 



3, 



to compute the proper- 
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221. The Mie fluid 



ties of the Mie(n,m) model fluids 9 
potential is often said to be short-ranged, however, for any finite system size, the potential is 
not rigorously zero at a distance of the half box length where the potential is typically trun- 
cated. Some approaches were proposed for dealing with the long-range tail of the potential 
0, Q. Perhaps, due to the commonly used procedure of potential cut-off during molecular 
simulations, there are lot of contradictions in the literature. As an example, the recent 
results of Okumura and Yonezawa [18( have indicated that the coexistence curve for the Mie 
(n,6) models scaled by the critical temperature and density practically coincides with each 
other, when 7 < n < 32. Dunikov and co-workers 9] have also demonstrated that the coexis- 
tence and interfacial properties of the LJ fluid approximately follows the corresponding state 
(CS) principle, if calculated with different cut-off's of the interaction potential. Meanwhile, 



Galliero and co-workers 



121 ]. by studying modified forms of the Mie(n,m) potential, found 



that no CS approach is possible between potentials having different repulsive exponents as 
well as different functional forms. 

It is well accepted that CS rule permits the prediction of unknown properties of many 
fluids from the known properties of a few [l, [2J, ! 
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261 ] . Its application to the model 



potential fluids allows to avoid the usually timeconsuming molecular simulations, and also 



makes possible the utilization of some important developments in statistical mechanics which 
would otherwise be prohibited by computational difficulties. Firstly the CS rule was derived 
by van der Waals based on his well-known equation of state, 



f( Pr ,T r ,Pr) =0, (2) 

where the variables were reduced by their critical values, p R = p/p c ', (reduced number den- 
sity or inverse molar volume), T R = T/T c (reduced temperature), and P R = P/P c (reduced 
pressure). Later, the CS principle was extended by introducing additional parameters to 
the study of various types of molecular fluids Thus, in general, macroscopic CS law 
states that all substances obey the same equation of state in terms of the reduced variables, 
or, in other words, the state of a system may be described by any two of the three variables: 
pressure, density, and temperature. 



Very recently [27[ , we have studied the application of the CS principle for the hard sphere 

2811 . neither 



attractive Yukawa (HAY) fluid. This model fluid is far from being conformal 
its reduced second virial coefficient B% evaluated at the critical temperature for various k 
(range of attractive tail) is constant. However, we have shown, that for different values of 
K, HAY fluid obeys the CS law for various thermodynamic properties. Unlike the studies 
of Refs. 
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28 



29| where the approximation of the constant value of B\ was used, we 



have proposed a new rescaling of fluid density, that allowed us to obtain a relation like 
eq.([2]) for the HAY fluid. In Ref. 27J we have suggested that the same kind of rescaling 
may be universal for different pair potentials, what afterwards has been partially confirmed 
17| . The purpose of the present investigation is to show, by applying canonical Monte 



Carlo (MC) simulation, that different thermodynamic properties of t 



re Mie fluid obey the 



principle of the corresponding states. Unlike the previous studies [12j, [18|, here we present 
the simulation data not only for the vapor-liquid coexistence densities, but also results of 
the surface tension, interfacial width, and supercritical pressure calculations. 
In this work we study the systems with two particular cases of the potential ([T]). Namely, 
the one-parameter Mie(2m,m) model, in which increasing m leads to a shorter ranged 



potential. Our calculations are focused on three values o 
studied here is Mie(n,6). In this case, like in Refs. [12l . 



m: 7, 9 and 12. Another potential 
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19l . l2l| , only the repulsive part 



is manipulated; three potentials with n = 12, 18, and 32 are considered. 



The applied simulation method is the same that was used in our previous works 27|, |30j, 131 ]. 
Briefly, the simulations of the vapor-liquid interfaces were performed in a parallelepiped 
cell with sides L x = L y = 12, and L z was chosen to be, at least, three times longer than L x ; 
the number of particles inside the box was iV > 1500, depending on the thermodynamic 
conditions considered. Periodic boundary conditions were applied in all three directions. 
The pressure of the supercritical Mie(n,m) fluid was calculated in a cubic simulation cell 
with L x = 12. We used r cut = 5.5 for Mie(12,6) and r cut = 5.0 for other systems. As 
we verified, these values of r cut are sufficiently long to omit the influence of the potential 
truncation on the averaged results presented below. 

Coexistence vapor p v and liquid pi densities, and the interfacial width, 5, were obtained 
at the end of each simulation run by fitting the density profiles, p{z), to the following 
hyperbolic tangent function 

1 1 /2(z — z ) \ 
P\ z ) = 2^ pL + pv > ~ 2 tanh { § J' ( 3 ) 



where z is the position of the Gibbs dividing surface. The critical density and 
for these model fluids were calculated by using the rectilinear diameters law 



temperature 



32| with the 



universal value of (3 = 0.325. The critical pressures were estimated on the base of the 



Clausius-Clapeyron equation [l]. Critical parameters are given in Table I. 

Fig. 1. depicts the reduced density, pn, as a function of the reduced temperature, Tr. It 
is clearly seen that the coexistence curves of Mie(2m,m) fluid with m = 7, 9 and 12, and 
Mie(n,6) fluid with n = 12, 18, and 32 map onto a single master curve. As was mentioned 

nn 

above, the similar behavior for Mie(n,6) and LJ fluids has been reported in Refs. [9j, 118] . 
These results confirm that the vapor-liquid phase diagrams of any Mie(n,m) potential obey 
the CS law, at least in the range of considered parameters. 



As usual 3jJ, the interfacial tension is calculated by 



7= ^{{P. Z ) ~\m,) + (Pyy)]}, 



(4) 
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where, Pa are the components of the pressure tensor. 

Our results of the surface tension are presented in its reduced form 



23|, 



* = - vL- (5) 

pc T c 

The values of 7 r as a function of the reduced temperature Tr are plotted in Fig. 2. As seen, 
the six sets of the surface tension data form a single master curve, which means that the 
surface tension of Mie fluid also obeys the CS theorem. In other words, application of the 
corresponding state rule may avoid the time consuming calculations of the surface tension 
for other pairs of parameters m and n if the critical parameters are known. It is worth 
noting that the master curve in Fig. 2 is slightly different from its HAY fluid counterpart 



reported in Ref. 



271 ]. and is represented by the following empirical equation, 



7r = 4.8(1 -Tr) 11 / 9 , (6) 

which resembles the expression proposed by Guggenheim to correlate the experimental data 
23]. 



In many engineering problems (for example, for dispersant applications in petroleum and 
pharmaceutical industries) it is of high importance to predict the vapor-liquid interfacial 
width as a function of temperature. Since such calculations are not trivial due to the spacial 
fluctuations of the interfaces, application of the CS principle might be of high utility. Since 
the value of surface tension is inversely proportional to the interface width [3, 0], it is 
natural to expect that the rescaling of 5 must be like Eq.Q, 

Pc T c 

In Fig. 3 the reduced widths of the vapor-liquid interface, S r , as a function of Tr are 
shown. As seen, the interface width curves of the Mie fluids are almost overlapped with 
each other, taking into account the error bars. It confirms that calculation of interfacial 
width can also be simplified by application of the CS principle with rescaling ([7j). Moreover, 
our preliminary calculations of 5 r for some parameters of HAY fluids indicate that it is 
practically the same as the one presented in Fig. 3. 
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Finally, we analyze the reduced pressure Pr as a function of reduced density p r , defined as 

Pr = ~§/3- ( 8 ) 
Pc 

Such rescaling of fluid density has been proposed and tentatively justified in our previous 
work [271 ]. where the extended corresponding states law has been applied for the description 



of the PVT properties the HAY fluid. 

The values of Pr for three reduced temperatures are presented in Fig. 4; as seen, the reduced 
pressure data for all Mie fluid potentials in question match on the same master curves at 
each value of Tr. Besides, the results of Pr for the Mie(n,m) family fluids and the HAY 
fluid [l7j practically coincide. The best agreement is reached at the lowest temperature and 
densities; only slight discrepancy, almost within the error bars, is observed at high fluid 
pressures, Pr > 45. The discrepancy is expected to be more pronounced at higher reduced 
pressures. By using the critical parameters presented in Table II, we also show in Fig. 4 the 
reduced pressure for the three real fluids, Ar,N 2 , and CO. As expected [15], the pressure 
of argon is predicted quite well by the Mie potential model, while some deviation of N 2 and 
CO pressures from the simulation data may be attributed to the slight non sphericity of 
these gas molecules. 

In summary, we have calculated the vapor-liquid equilibrium and interfacial properties of 
Mie(n,m) fluids using canonical Monte Carlo simulations. Our results indicate that the 
coexistence densities and the surface tension, as well as the width of the vapor-liquid interface 
obey the corresponding states principle. On the base of new accurate critical parameters 
(evaluated from our MC simulations), it is found that the critical compressibility factor is 



around 0.3, which agrees well with experimental data for some real substances 23j, |34|. The 
new density rescaling, p r = p/pl^ 3 , is universal for the calculation of the reduced pressure of 
Mie(n,m), HAY j^], and Sutherland fluids We expect that the proposed new criteria 
could be of relevant help for testing new theoretical approaches for the investigation of the 
del, as well as real fluid systems [3, [si, Q, 11. lo1|. 
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Figure Captions 

Fig. 1 Reduced vapor-liquid coexistence curves of the Mie(n,m) model fluid at different 
interaction ranges. The solid line is just a guide for the eye, which is well described by the 
empirical formulae proposed by Guggenheim [23J to build the coexistence curve for argon. 
Error bars do not exceed the symbol size. 



Fig. 2 Reduced surface tension 7 r as a function of reduced temperature Tr for the same 
systems considered in Fig.l. 

Fig. 3 Reduced interfacial width, 5 r , as a function of Tr for the same systems considered 
in Fig.l. The solid line is a fitting curve, 5 r = ^jo.se ■ 



Fig. 4 Reduced pressure, Pr, as function of reduced density, p r , at three reduced 
temperatures Tr = 1.23,1.52, 1.92, and 4.651 (from bottom to top). Our data for Mie 
(12,6) at Tr = 4.651 coincide with the results of Johnson et al. [2J. Symbols and lines 



depict the simulation data of Mie fluid, and experimental estimations for real gases[33j, |34j, 
respectively. Solid line is a guide for the eye. 
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TABLE I: Critical parameters of Mie fluids 



Systems T c p c P c Z c 



Mie(24-12) 


0.575 6 


0.3909 


O.O663 


0.294 7 


Mie(18-9) 


0.730 8 


0.360 9 


0.078 2 


0.297 8 


Mie(14-7) 


0.990 7 


0.330 8 


0.098 2 


O.3OO7 


Mie (12-6) 


1.2909 


0.314 6 


O.II83 


0.291 8 


Mie(18-6) 


1.050b 


0.3305 


0.102 2 


0.294 6 


Mie(32-6) 


0.865 5 


0.3407 


0.090 3 


O.3O65 



TABLE II: Critical parameters [331 ] and diameters [3j] a of the three fluids analyzed in Fig. 4 



Substances 


T C (K) 


p c (mol/L) 


Pc (MPa) 


cx(pm) 


N 2 


126.19 


11.1839 


3.3958 


375 


Ar 


150.69 


13.4074 


4.863 


342 


CO 


132.86 


10.85 


3.4935 


371 
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